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Abstract 

In the Maxwell-Chern-Simons theory coupled to Nf flavors of 4-component 
fermions (or even number of 2-component fermions) we construct the gauge- 
covariant effective potential written in terms of two order parameters which are 
able to probe the breakdown of chiral symmetry and parity. In the absence of 
the bare Chern-Simons term, we show that the chiral symmetry is spontaneously 
broken for fermion flavors Nf below a certain finite critical number Nj, while the 
parity is not broken spontaneously. This chiral phase transition is of the second or- 
der. In the presence of the bare Chern-Simons term, on the other hand, the chiral 
phase transition associated with the spontaneous breaking of chiral symmetry is 
shown to continue to exist, although the parity is explicitly broken. However it is 
shown that the existence of the bare Chern-Simons term changes the order of the 
chiral transition into the first order, no matter how small the bare Chern-Simons 
coefficient may be. This gauge-invariant result is consistent with that recently 
obtained by the Schwinger-Dyson equation in the non-local gauge. 
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1 Introduction 



The Chern-Simons term |jl|, ^ ^ generates a gauge-invariant mass for the gauge 
field, combined with the ordinary Maxwell term. Formally, in the Chern-Simons 
theory which does not include the kinetic term for the gauge field, the gauge 
sector of the theory is strictly renormalizable, although the Maxwell theory is su- 
perrenormalizable in 2+1 space-time dimensions. The Chern-Simons term affects 
the long-distance behavior of the theory and the perturbative expansion of Chern- 
Simons theory contains logarithmic divergences ^ 0. However, the topological 
nature of the theory should allow only trivial, finite renormalization [Q. 

Actually it is verified that this happens in the various renormalization scheme. 
The massive matter fields may have a finite one-loop radiative correction to the 
statistical parameter 9, the coefficient of the Chern-Simons term. It is generally 
known as the Coleman-Hill nonrenormalization theorem |^ that the topological 
mass term at one-loop order is the only one radiative correction and there are 
no radiative corrections from higher-order loops to the topological mass term. 
Therefore the coefficient of the topological mass term at one-loop order is exact 
to all higher orders. The requirement of this nonrenormalization theorem is the 
gauge invariance and analyticity of the one-loop n-point functions. 

Indeed, from the explicit consideration of the Feynman graph, the following 
results are obtained. When there is no topological mass at the tree level, no 
topological mass terms are induced radiatively to the two-loop order in either 
Abelian or non-Abelian gauge theories with massive fermions in 2+1 space-time 
dimensions P|. Moreover, if there is a nonzero tree-level topological mass, the 
radiative correction to the topological mass is shown to vanish in Abelian gauge 
theory coupled to the massive fermion 0. However it should be remarked that 
the existence of massless charged particles or spontaneous gauge symmetry break- 
ing violates the initial assumptions of the Coleman-Hill theorem and higher loop 
renormalizations of the Chern-Simons term appear | ]IU| , |TT], 0. 

The nonrenormalization theorem in the weak coupling expansion can be ex- 
tended to the large fiavor (Nf) case, i.e., to the 1/iV/ expansion: there is no 



radiative corrections to the Chern-Simons coefficient 6 in the 1/N expansion |T3 



This argument do not apply to the case in which the fermions remain massless 
due to absence of the analyticity. In the Abelian case, the nonrenormalization 



theorem to all orders is also shown by using the effective potential [0. It will be 
interesting to extend these analyses to the non-Abelian case such as QCD3 



le, 17 



In this paper we adopt the 1/A^/ expansion as a systematic truncation scheme. 
The 1/A^/ expansion is quite interesting from the viewpoint of the non-perturbative 
study, since a class of theories with a four-fermion interaction is shown to be renor- 



malizable in the framework of 1/A^/ expansion |jT8| in spite of its nonrenormaliz- 



ability in the weak coupling expansion. The Chern-Simons term may affect the 
high-energy behavior of the gauged four-fermion model [|1^, |19| and the pattern of 
dynamical generation for the fermion mass, although the four-fermion interaction 
is not investigated in this paper. 
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In this paper we regard the dimensional couphng constant e, the electric charge, 
as a genuine dimensional parameter that sets the natural scale of the theory. 
This viewpoint has been taken by Appelquist et al. PDI and several authors. 



The dynamically generated fermion will have a mass which is proportional to 
this scale: m^^ ~ e^. The fermion mass is generated when > 1/^/ where 

N'j is called the critical number of fermion flavors. Another way of viewing the 
dimensional parameter is to take it as the ultraviolet (UV) cutoff of the theory 
where A is the UV cutoff and g = e^/A is a dimensionless coupling [|21|. In 



this case, the only dimensional parameter in the model is the UV cutoff. In 
this standpoint, if the dynamical mass is generated, it will be independent of the 
UV cutoff A. This leads to the non-perturbative renormalization for g, i.e., the 
coupling is running, i.e., (3{g) := A^ 7^ 0. Therefore, the first viewpoint is not 
compatible with the second one and two viewpoints are in many ways different 
from each other. 

The Maxwell- Chern- Simons theory ^ |^, ^, ^ ^ with Nf flavors of Dirac 
fermions in (2+l)-dimensions is defined by the following euclidean Lagrangian 
density: 

^ = ^^Iv + -;^^iiup^tiduAp + Cgf 

+^°(z7;.5^ - me - m^r)^" + c^^Tm^^^m. (1-1) 

where the index a denotes the fermion flavors (a = 1, Nf) and Cgf is the gauge- 
fixing term specified in the next section. In the above Lagrangian density we have 
introduced an interpolating parameter (3 between the Maxwell (/3 = 1, ^ = 0) and 
the Chern-Simons (/3 = 0, 6' 7^ 0) theories. Our analysis based on the effective 
potential in this paper can be applied equally to the case of /5 = as well as /5 7^ 0. 
The (2+l)-dimensional QED (QED3) with a Chern-Simons term is a special case 
of this model. For the Chern-Simons theory, the term f -F^j, can be interpreted as 
a regularization which is removed finally, i.e. /3 — >■ in this case. 

We start from the 4-component formulation for fermions in (2+l)-dimensions 
[ pO|] . Each 4-component fermion [| \1/ is decomposed into two 2-component fermions 
il)i,il)2 as 

Then the Dirac conjugate \1/ := \l/^7o is written as 

^ = (^1,-^2) (1.3) 

by using ipa '■= ipla^la = 1, 2), see Appendix A. Accordingly the gamma matrices 
satisfying the Clifford algebra {7^,7;^} = — 2(5^j, are written by 4 x 4 matrices. In 
this reducible representation, there exist 4x4 matrices anticommuting with all 



^ We use 4* without the index a to denote each 4-component fermion. 

^In what foUows we use the notation ^4 _B in the sense that A is defined by B, while A = B 
impUes that A is identicaUy equal to B. 
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the gamma matrices, say 73,75, in sharp contrast to the irreducible representa- 
tion corresponding to 2-component fermions |20|. Then we can define the chiral 
symmetry similarly as in the (3-|-l)-dimensional case: 



\1/ exp(i7973)\l>, \1> \l>exp (2^973), (1.4) 

under which 'l'27^(9^\l/ is invariant. Moreover, we can define the parity transfor- 
mation: 

^'(xo,x) ^ V^ixo,xp), 

A^{xo,x) {-lY^'A^{xo,xp), (1.5) 

where xp = {—Xi,X2) for the spacial coordinate ) and V = -i757i- 

The parity transformation ( p..5|) is equivalent to ipii^O:^) ~^ O'i4'2{xo,xp) and 
i/j2{xo,x) aiipi{xo,Xp), for 2-component fermions. The transformation prop- 
erty of each term in the Lagrangian under two transformations is summarized as 
follows. 



Transformation property 


term 


chiral 


parity 




odd 


even 




even 


odd 


f: A A 


even 


odd 



The order parameter for chiral symmetry breaking is defined by 

(^r^) = {tPM - (1.6) 

while the order parameter for parity breaking is defined by 

(M) = (^1^1) + (V'2V'2). (1.7) 

According to analytical and numerical studies |Q of the Schwinger-Dyson 
equation for QED3 with a Chern-Simons term {j3 = 1,6 0,me = rrio = 0), we 
can write the schematic phase diagram for the Maxwell-Chern-Simons theory with 
Nf fiavors of 4-component fermions, see Fig. 1. First note that, when 6' 7^ 0, the 
parity is explicitly broken, {^t'^) 7^ 0. The solid line in Fig. 1 denotes the critical 
line of chiral phase transition below which the chiral symmetry is spontaneously 
broken, ('I'\E') 7^ 0. In Fig. 1 we have chosen typical five points, assuming that the 
critical number of fiavors is nonzero and finite: < Nj < 00. The points A 
and E in Fig. 1 are in the chiral symmetric phase; the point A is in the parity even 
phase, while E is in the parity odd phase. The point D is located just on the critical 
line. The points B and C are in the chiral-symmetry-breaking phase; the point 
B is in the parity even phase, while C is in the parity odd phase. For the order 
of phase transition, it is well known that the point {Nf, 9) = {Nj, 0) exhibits the 
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continuous second order transition through the analyses of the Schwinger-Dyson 



equation . 

By using the Schwinger-Dyson equation, the critical number of flavors N'j has 
been evaluated. Finite and nonzero critical number, < A^^ < oo has been 
reported: N'j = ^ = 3.2 first in the Landau gauge by Appelquist, Nash and 
Wijewardhana p^, and iV? = ^ = 4.3 in the heuristic treatment of the non-local 



/ ~ 37r2 

gauge by Nakatani and N'j = 3 ~ 4 by Nash |^ . These results are confirmed 
by the systematic investigation of the non-local gauge by Kondo et al. [^] and by 
Kondo and Maris On the other hand, infinity of Nj, Nf = oo, was claimed 
by and subsequent works [^9|, |3^ . Different approximations to the vertex 



function and the existence of infrared cutoff may lead to different results, see |3T 
and |S^. Such a controversy exists |Q ^ also in the approach of the Cornwall- 
Jackiw-Tomboulis (CJT) effective potential [^. The Monte Carlo simulation of 
non-compact lattice QED3 shows Nj = 3.5 ± 0.5 according to Dagotto, Kogut 
and Kocic [^]. The result are controversial up to now. This shows that it is 
extremely difficult to find a consistent truncation scheme in the Schwinger-Dyson 
approach which is able to deduce the gauge-invariant result. For more details on 
the Schwinger-Dyson equation of QED3, see for example a review ||37| . 

In the absence of the bare Chern-Simons term, the fermion mass is dynamically 
generated for Nf < N'j in such a way that the parity is preserved, i.e., (f'\E') ^ 
and (^'r\l/) = 0. This is rephrased in terms of the fermion dynamical mass mi, m2 
for two 2-component fermions, respectively. Here the chiral symmetry 

breaking, (^'\E') ^ 0, implies mi ^ m2 and no parity breaking, {^t"^) = 0, 
implies more specific pattern for the fermion dynamical mass: 



■m,m, ...,m, — m, — m, —m , 

Nf Nf 



This result is consistent with various analyses so far |38|, |41|, ^ 

as the general consideration by Vafa and Witten |]45 |. 
A new and remarkable feature discovered in 



as well 



^ is that the order of chiral 
phase transition turns into the first order when 6^0. Namely the chiral phase 
transition is of the first order on the whole critical line except for one point 
{Nf, 9) = {Nf, 0) in the phase diagram, which is in sharp contrast to the previous 



analysis The meaning of the first order transition is explained from the 

behavior of dynamically generated fermion mass as shown in Fig. 2. In Fig. 2, 
fermion masses mi,m2 are shown as different functions of 6 when Nf < Nf. 
At the point B, 6 = and fermion masses for two 2-component fermions are 
generated in such a way that the parity is preserved: mi = — m2 > (we choose a 
convention such that mi > without loss of generality) and the chiral symmetry 
is spontaneously broken, (^'^) ^ 0. When 6 ^ 0, we see this pattern for mass 
generation is destroyed so that |mi| 7^ |m2| (mi > 0,m2 < 0). As 9 is increased 
while Nf is kept fixed, each ma{a = 1,2) increases monotonically. However, at 
9 = 9c, m2{9) jumps discontinuously towards mi{9). For 9 > 9c, it is realized that 
mi{9) = m2{9), which implies restoration of the chiral symmetry: ('I'\E') = 0. In 
this sense the chiral phase transition is of the first order when 9^0. 



5 



However this result is available only in a specific gauge. In the Schwinger- 
Dyson equation approach ||2^, the non-local gauge was adopted in order to avoid 



the vertex correction. A possible way to improve the gauge-dependence in the 
Schwinger-Dyson equation approach is to improve the vertex so as to be consistent 
with the Ward-Takahashi identity. However this way is not unique. 

The purpose of this paper is show the existence of the first order phase tran- 
sition in the presence of the bare Chern-Simons term as well as the second order 
transition in the absence of the bare Chern-Simons term in the gauge-invariant 
way. For this, we construct the gauge- covariant (gauge-parameter independent) 
effective potential in terms of two order parameters: = {^'^) and x = (^^t\1/), 
up to the leading order of expansion. To construct the effective potential, 
we adopt the inversion method PSI . 



This paper is organized as follows. In section 2, we show how the gauge- 
covariant result is obtained for the gauge-invariant order parameters through the 
evaluation of the effective potential in the scheme of l/Nj expansion. In section 
3, after introducing two external sources which trigger the spontaneous breaking 
of the chiral symmetry and parity, we show how explicitly the order parameter 
is related to the vacuum polarization of the photon propagator. In section 4, we 
observe the dependence of the order parameter on two external sources in detail. 
Sections 2 to 4 are preparations for performing the inversion. Now, in section 
5, we state a general strategy of obtaining the non-trivial solutions through the 
inversion method in the presence of two external sources. In section 6, as a 
special case we study the system with a single order parameter and in a slightly 
simplified situation we show that the chiral symmetry is spontaneously broken in 
the three-dimensional gauge theory, e.g., QED3 by using the inversion method. 
In section 7, we give the argument on the behavior of the effective potential near 
the stationary point in the absence of the Chern-Simons term. The stationary 
point corresponds to the non-trivial solution (order parameter) corresponding to 
the spontaneous breakdown of the relevant symmetry. In section 8, we study the 
effect of the Chern-Simons term on the chiral phase transition by extending the 
effective potential given in section 7. The final section is devoted to conclusion 
and discussion. In Appendix A, we give a formulation of 4-component fermion 
in (2-1-1 )-dimensions. In Appendix B, we derive the gauge covariant formula in 
the 1/Nf expansion given in section 2 for the gauge-invariant order parameter. In 
Appendix C, the vacuum polarization tensor in three-dimensional gauge theory is 
obtained together with its power-series expansion in the external source. 

2 Order parameter in the l/Nf expansion 

In this paper we take the covariant gauge-fixing: 

Cgf = ^{d,A,r. (2.1) 
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From the Lagrangian ( |1.1| ) with the gauge- fixing term ( |2.1| ), it is easy to see that 
the inverse of the bare photon propagator reads 



(2.2) 



and hence the bare photon propagator is given by 



+ 02 



5. 



k^kjy 



+ 



k/^ki/ ^ (^/lupkp 
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(2.3) 



where we have used eoi2 = 1 and e^^x<^pax = 6^p6^^ - 6^^6up. 

The general form of the vacuum polarization tensor in three-dimensional gauge 
theory is written as 



Ilpuik) 



k^ky 



IiT{k) + ep,p^Iio{k). 



(2.4) 



The full photon propagator Dp^{k) is related to the vacuum polarization tensor 
Il^^{k) and the bare photon propagator Dj^J{k) through the Schwinger-Dyson 
equation: 



Dp^'\k) = D<^t\k)-UpAk). 
Therefore the full photon propagator is written in the form 

D^uik) = 



A;2 



DT{k) + ^D,{k) + e,^p^Do{k), 



(2.5) 



(2.6) 



with 



Drik) 
Do{k) 



(3e - Urik) 



[(3k^ - nT(fc)]2 + [Qy/k^ - Yio{k)f ' 

eV¥-Uo{k) 



[f3k^ - Urik)]^ + [eVk^ - Hoik)] 



(2.7) 



In order to calculate the order parameter, we introduce the source term 

Cj = Jefr"(a;)^"(x) + Jof'"(x)r*"(a;). (2.8) 
Then the partition function in the presence of the source is defined by 



Z[Je, Jo] = J V'^V'^VAf.exp - J d^x{C + 

and the generating functional is given by 

W[J,,Jo\ = -\nZ[J,,Jo]. 



(2.9) 



(2.10) 
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The order parameter is obtained by differentiating the generating functional with 
respect to the corresponding source. For the translation-invariant order parame- 



ter, we obtain 



1__5_ 



W[J,, Jo] 



Je = Jo=0 



Je = Jo = 0' 



(2.11) 



where ^ = J d^x is the space-time volume. The generating functional is given by 
the sum of all the bubble diagrams which is depicted in Fig. 3 up to the leading 
order of 1/A^/ expansion. Separating the bare (interaction free) part with the 
superscript (0), the order parameter which we calculate can be written in the 
form (see Appendix B): 



(0) 



(27r)3 



Here n|^^^ denotes the vacuum polarization in which all the radiative corrections 
up to m-(fermion) loops are included, and D^jf) denotes the photon propaga- 
tor which is constructed from Eq.( |2.5|) by using H^"^. We use the convention: 



uPik) 



= IlQ\k) and D^^J is the bare photon propagator. 



To carry out the calculation in the 1/A^/ scheme, therefore, we need to calculate 



d^k 



(2vr) 



(27r)3 \ai.2 _ T^(■^)^u\^2 i \q./u2 _ rrW^ 



- U^ik)]^ + [eVk^ - W^o'ik)] 



(2.13) 



where we have used e^^^^^^is = 26ai3- 

Up to the leading order of 1/Nf expansion, we have only to consider the case: 
n = l,m = 1, which is shown in Appendix B. This is a novel observation in 
this paper. This should be compared with the usual perturbation theory with 
respect to the coupling constant e^. In this case, the order of the graph is counted 
with respect to the coupling, e^. In the lowest order of the perturbation series 
with respect to the coupling constant, we need to choose n = 0,m = 1 as is 

El- 



It should be remarked that 



well known in the four- dimensional case D = A 
the expression (|2.13| ) is free from the gauge-fixing parameter ^. Therefore, if we 
calculate the gauge-invariant order parameter according to this expression, we will 
get the gauge-covariant (i.e. gauge-parameter independent) result. 



8 



This should be compared with the Schwinger-Dyson equation approach in 
which a specific gauge has to be chosen even in calculating the gauge-invariant 
quantity. The Schwinger-Dyson equation approach does not in general guarantee 
that the gauge-invariant quantity calculated in one gauge coincides with that 
calculated in other gauges. 

Even in the approach of the effective potential, a special gauge is usually taken 
such as the Landau gauge, see |4^. Note that if we are allowed to change the 
order of calculating the graph or equivalently performing the integration over the 
internal momenta, we get (see Fig. 3(b)) 



/ (0^- = / |^Sl»>(p)E<')(p). (2.14) 



In this case, even if we start from the second expression, we should obtain the 
same result as long as it is finite. However, it is divergent and hence we need to 
regularize it. The naive cutoff procedure usually taken in the calculation of the 
effective potential may break the gauge invariance in the sense that this reordering 
is not allowed and hence the result may be not gauge- covariant. Therefore, if we 
start from the gauge- noninvariant quantity S(p), it is quite difficult to recover the 
gauge covariance after the calculation. 

3 Order parameter and vacuum polarization 

We introduce the projection operators 

X±:=\{l + r\ r=(^ _° J , (3.1) 

which satisfy the properties: x± = Xi? X+X- = ^ind x+ + X- = 1- Then the 
source term (E.8|) is rewritten as 



Cj = J+^"(x)x+^"(x) + J_^"(a;)x-^"(x), (3.2) 
where we have defined 



or equivalently 



J± := Je ± Jo, (3.3) 



Je = \{J+ + J-), Jo = \{J+-J-). (3.4) 



3.1 free part 

The bare fermion propagator S^'^'^ (p) is also decomposed as 



^Note that the signature of J^, is opposite to that of the bare fermion mass mg, TOo. 



9 



where 



1 



Then the free parts of two order parameters are decomposed as 
with 



(0) 



By introducing the UV cutoff Af, we have 



(0) 

± 



/ dp p J± 

TC"^ p'^ + Jj.' 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



Here it should be remarked that the quantity defined by eq. (|3.9|) has 

an ambiguity arising from the lower bound of the integration, no infrared cutoff, 
unless we specify the signature of J±. In order to define this quantity unambigu- 
ously, we introduce the small and positive infrared cutoff e > and then take the 
limit e I 0. Thus we obtain 



(0) 



lim 
lim 



dp- 



p'J± 



p^ + Jl 
J^{Af-e) + J[ 



J± J± 
arctan arctan ■ 



A/ 



J±Af + ( arctan — —sgn{J±) 



J-. 
A/ 



71 

2' 



where sgn{J±) denotes the signature of J±. Therefore we have 



(0) 



A2 



A/ 



(3.10) 



(3.11) 



In what follows, we assume that J+ is always positive, J+ > 0, without loss 
of generality. Then, depending on the signature of J_, there are two possibilities: 
Je > Jo for J_ > and Jg < Jo for J_ < 0. Here it should be remarked that the 
signature of Jg and Jo are not specified. 

First we consider the case (I) in which J+, J- > 0, i.e. Je > Jo- Thus we have 



(0) 



A2 



(0) 



2 
1 
2 



J+ + J- TT Jl 

A/ " 2~ 

Je VT J2 + J2 



J2 



A2 



0(J^ 



A 
"J4 



/ 



2 



A2 



TT 



C(J3) 

J2 



A 



Jo 

A/ 



TT 2JeJo 

2^ 



A2 
+ 0(J3 



+ C(J' 



(3.12) 
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In the second case (II) where J+ > 0, J_ < 0, i.e. Jo > Je- we obtain 



(0) 



2 



A/ 



2 



A 



/ 



2 A 



J, - J_ 



A, 



/ 

2 A^ 



Jo 

A/ 



TT 72 I 72 



(3.13) 



3.2 vacuum polarization 

By using the decomposition ( p.5| ) for S^'^\p), the vacuum polarization tensor at 
one-loop is also decomposed as 



(27r)3 

iV^[n«(A:;J+) + n«(A;;J_)] 



(3.14) 



where we have defined 



(2vr) 



For the one-loop vacuum polarization tensor with the tensor structure: 



nS(fc) 



kr, 



(3.16) 



we find 



U^^\k) = Nf[U^fi{k;J+)+U^^\k;J.)], 
ng)(fc) = Ar;[ng)(A;;J+)-ng)(A;;J.)], 



(3.17) 



since the tensor Il^^}{k; J±) should have the same tensor structure as Il^^^{k). Note 

that there exists a minus sign in IlQ\k). Actual calculation shows (see Appendix 
C) that 



n5!^(fc; J) 



in 



■ , fc2 - 4J2 k 

2 J H : arctan 



k 



2\J\ 



and 



(1) ^ k 

Yin {k\ J) = — J arctan— — - 
o ^ ' ' 27r 2 J 



(3.18) 



(3.19) 
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It should be remarked that Il^j^\k; J) is an even function of J and IlQ\k; J) is 
odd in J. Therefore, if Jg = 0, then J± = ±Jo and 

U^T\k) = 2NfU^T\k; Jo), Tl^oik) = 2NfU'o\k; Jo). (3.20) 

On the other hand, if Jo = 0, then = = and 

U^^\k) = 2NfU^T\k; Je), Il'b\k) = 0. (3.21) 

The usual criterion for generation of the induced Chern-Simons term is given 
by Ilo\k 0) 7^ 0- The induced Chern-Simons term is generated when > 
0, J_ < 0(Je < Jo), since 

ng^(fc-O) = ^[sgn{J^) - sgn{J_)] 



An 

> 

-[sgn{Je + Jo) - sgn{Je- Jo)] 

4:71 

JO (Je> J.;J+,J->0) 

- (Je< J.;J+>0,J_<0) • ^'-''^ 

4 Order parameter as a function of source 

We introduce new mass parameters (or sources) Ji, J2 corresponding to two 2- 
component fermions, ipi^ip^'- 

Ji = J+, J2 = -J„, (4.1) 
so that the source term ( p.8|) is rewritten as 

Cj = JiV'i(x)^i(x) + J2^2{x)i^2{x), (4.2) 

where 

Je=\{Jl-J2), Jo = \{Jl + J2), (4.3) 

or 

Jl = Je + Jo, J2 = - Je + Jo. (4.4) 

Similarly in the preceding section, the free parts of the order parameter are 
decomposed as 

(^^^)(0) ^ (^i^i)W + (V^2^2)^°\ (4.5) 



where we have defined: 



{i'ai^ar= t'%4^2 («=1'2). (4.6) 
JO TX'^ + 
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We apply the above decomposition also to the interaction part. By taking into 
account eq. ( p.l2|) and 



d 



d 



d 



d 



d d 

+ 



dJe dJi (9 J2 dJo dJi dJ2 
the two order parameters are rewritten in the form: 



(4.7) 



(4i 



where 



-1 



d 



i2n) 

X { - n? (k)] + [9V¥- n« (A:)] (k) 



dJa 



dJa 



+0 {{i/Nfy 



(4.9) 



Now we study the dependence of the integrand of Eq. (|4.9| ) on the source. 
First of all, it is shown by the same procedure as in the previous section that the 
free part ( [4.6|) reads 



(0) 



Ja 



TX 



J' 



(4.10) 



where cTq denotes the signature of Ja, i.e., ctq = jj^. 

To study the dependence of the interaction part on the source, we observe that 
the one-loop calculation in Appendix C leads to 



U^^\k) 

ng)(fc) 



ak 
2a 



P '^9 rr 73 

V fc2 377 V fc3 



n n k 



+ o (.Jt) 



(4.11) 



and 



any^(fc) 

dJa 

dJa 



4a 



Ja 
k 



7T k"^ 



2a 



8Ja 

'tt k 



+ o{Jt) 
+ o{j!), 



where a parameter a is defined by 



Nfe"^ 



a := 



(4.12) 



(4.13) 



which has the dimension of mass in three- dimensions. 
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4.1 case (I) 

First we consider the case of Ji > 0, J2 < (J+ > 0, J_ > 0). The numerator of 
the integrand in eq. ( [4.9|) has the following dependence on the source J a'- 



where 



Po{k) = 2eak, 
16 

P^{k) = -aa—9a + 4a{pk + a), 

TT 

p^(k) = -W < 0, 

TT K 

00 2 
TT k 

On the other hand, the denominator of the integrand in eq. 

{[(3k^ - n«(A;)]2 + [6V¥ - n«(A;)]2}^i 



(4.15) 



) has 



[F{k)]-' 



+ 0{J!), (4.16) 



where 



Qoik) 
Qi{k) 
Q2{k) 



iOak 

m 

-Aa 
4a 



40 



F{ky 

Pk + a 



-Aa 



^[{(3k^ + akf + {ekf]-Ae^k^ 



F{ky 



(4.17) 



and 



F{k) := + akY + {ekf. 



(4.18) 
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Substituting eq. ( [4.14| ) and eq. ( ^.16| ) into eq. ( |4.9| ), we obtain for a = 1,2 

b b b b . 

+0 , (4.19) 



where 



A{k) = Po{k)Qo{k) + P2{k), 
B{k) = Po{k)Qi{k) + P,{k), 
C\k) = PmQo{k) + P^{k), 
D{k) = Po{k)Q2{k), 

E{k) = Poik)Q,{k) + PMk). (4.20) 
Thus we obtain the order parameter in the form: 



^1 
where 



:= (^r^n/A^/ = Ve + u^Ji + vJ2 + Si(Ji)2 + hiJ^f + ri Ji J2 + O (j^) , 
:= (^2>2 = Ve + vJi + U2J2 + Wif + S2{J2f + r2 Ji J2 + O (|l)2^) 



Ua = 



1 r d^k 

sa = j^j ^[^3 (^) + ^aB{k) + C"^(fc) + D{k) + i?(A:)]/F(A;), 

ta = J^J -^,[-CJaB{k)+D{k)+E{k)]/F{k), 

ra = — ^^[C'^{k) + 2E{k)]/F{k), (4.22) 



and 



It is easy to check that 

n = 2t2, r2 = 2h, (4.24) 
since P5 = 2P4 and P{^go = 2Po(Q2 + cXaQi). 
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The parity even and odd order parameters are written as 



:= (*V^'")/iV/ = ¥^+-¥^- = (^i + (^2. (4.25) 

They obey 

(/7e = {ui^U2-'2v)Je^{ui-U2)Jo 

+ - S2 + 3(^1 - t2)]Jl + [Sl - S2 - (^1 - t2)]Jl 
+2[si + S2 - (tl + t2)] JeJo + C)(J^), 

(/?o = 2(/?e + (lii - U2)Je + (Mi + M2 + 2^)70 

+ + S2 - {ti + ts)] + [Sl + S2 + 3(ti + ^2)] Jo 

+2[si -S2- (h - t2)]JeJo + 0{J^). (4.26) 

In order to study the small 9 case, we first observe that 

A{k) = P2{k) + 0{e^) <Q, 
B{k) = P4{k) + o{e^) > 0, 

m = -^^ + 0(^^), (4.27) 



where 



Fo{k) := (/3A;2 + akf. (4.28) 



Note that when ^ = 0, Pq = 0, Qo = , Qi = and P^ = Therefore, the 
functions Qai^ — 1,2,3) do not contribute to the order parameter. 
Thus the coefficient is calculated up to 0{1/Nf) and 0{9): 



dp 1 f d^k Aa(3k 



Sa 



7r2 NfJ (27r)3Po(A;) 

1 /• d^A; (7jP3 + P4(A;) + P5(A;)] 



TV^y (27r)3 Fo(A:) 

^ r d'k 2Wpk r 
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■fj (27r)3 F,{k) 



+ 



f d^k Sa^k{/3k-2a) .^(n2,.r2\ oq^ 



and 



d^k 2ak 



NfJ (27r)3Fo(A;) 
It is easy to see that 

U1-U2 = 0(e/Nf), ui + U2^(ui + U2)e=o + 0(eyN]), 

s, + s2 - 0{e/Nf), si-S2^(si-S2)e=o + 0{eyN]), 

h + t2 = o{e/Nf), t,-t2 = {h-t2)0=o + o(eyNj). 



(4.30) 



(4.31) 



Hence, for example, we obtain Ui + U2±2v — 2(u0=o ± vo=o) + O (O'^/Nj^ . 

Finally, we consider the limit ^ — 0, i.e., the Chern-Simons term is absent. In 
this limit, we find Si + S2 — * and ti + 12 ~* 0, and Ui, U2 approach to the same 
limit: Ui,U2 u. Thus, in the limit 6' 0, two order parameters obey 

(^e = 2{U-V)J,+ [S^- S2 + 'i{h-t2)]Jl + [si- S2-{h-t2)]Jl + 0{J^), 

ifo = 2{u + v)Jo + 2[s^-S2-{h-t2)\JeJo + 0{j'^), (4.32) 
where 



u 



tn 



_ r^f dp 1 r d^k Aa(3k / 1 \ 
h ^^WfJ (27r)3 Fo{k) ^ \N])' 

- -I 
NfJ 

NfJ 



d^k 4q;2 / 1 ^ 

Hni 



{2nyFo{k) 

d^k a a 



16a 




TT 


h 


16\ 




TT / 


k 



(4.33) 



For example, if we introduce another cutoff Ap, we get 

2A/ ' 



2{u - v) 



(4.34) 



where 



~ A/ 



TT^ /-Ap d^k 4a{/3k + a) 



(27r)3 {(3k^ + aky 



a 



In 1 + 



(4.35) 
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4.2 case (II) 

Next we consider the case of Ji > 0, J2 > (7+ > 0, J_ < 0). 

In this case we take cTq = 1 for a = 1, 2 in Eq. (|4.11|) and Eq. ( [4.12|) . Then we 
obtain 

(^1 = ipe + u.h+vJ2 + s{.hf + t{J2f + 2tJiJ2 + 0(^.fi) , 

^2 = ipe + vJi + uJ2 + t{Jif + s{J2y + 2tJ,J2 + 0(j^) , (4.36) 



where the coefficients u,v,s,t are given by Eq. ([4.29|) with 0"^ = 1. Hence the 
order parameter is written in the form: 



= 2{u-v)J, + 4{s-t)JeJo + 0{J^), 
ifo = 2ipe + 2{u + v)Jo + 2{s-t)J^ + 2{s + 3t)J^ + OiJ^) 



(4.37) 



When 6' = 0, for example, we get 



2(u + v) 



2k t 



(4.38) 



where 



2 d^k 4a{f3k-a) 
f ' 



(27r)3 (/3fc2 + aky 

-1 



A/ 



a 



In 1 + 



PA, 



a 



/3Ap 



1 + 



/3Ap 



(4.39) 



5 Inversion method 

In this section we restrict our consideration to the case in which the bare Chern- 
Simons term is absent. The effect of the bare Chern-Simons term will be inves- 
tigated in section 8. In this paper we use the inversion method to obtain the 
effective potential and find the non-trivial solution as its stationary point. The 
inversion method was proposed by Fukuda |4^ and was applied to gauge theories: 



for example, QCD4 [^], the strong coupling QED4 the gauged Yukawa 



model ^ and the Thirring model (reformulated as a gauge theory) |Q 



5.1 general strategy 

Let g denote an expansion parameter in a certain scheme of the expansion. For 
example, g may be the coupling constant in the usual perturbation theory 
or the expansion parameter in the expansion. In order to calcu- 

late the order parameter ('I'\E') and/or ('I'r\E') we introduce the source term 
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Lj = Je^"(a^)^"(a;) + Jo^E'"(x)r\l/"(x) into the original Lagrangian. We define 
the dimensionless order parameter and x by 

^ 2A} Nf ' 



vr2 (^"r^") 

with an appropriate factor which simphfies the following calculations. In a certain 
expansion scheme, suppose that the order parameter is calculated in the power 
series of g: 

oo 

= he[Je,Jo\ = ^g"4>n{Je,Jo), 
n=0 
oo 

X = K[Je,Jo\ = T.9''XniJe,Jo). (5.2) 

n=0 

Then we invert these equations with respect to the source: 

oo 

Je = /e[0,X]= E^?"^/m(0,X), 
m=0 
oo 

Jo = fo[<P,x]=T.9'^fMx)- (5.3) 

m=0 

We call this procedure the inversion. Here the coefficient functions fmi'PiX) ^i-nd 
fmi'Pi x) ci^re determined from the requirement that the following consistency con- 
dition should hold: 

= ^e[/e[0,X],/o[0,X]], 

X = ho[fe[(P,X]Jo[4>,X]], (5.4) 

and are written by using (pn{Je, Jo) and Xn{Je, Jo)- In the inverted series eq. (p. 3D, 
we look for the non-trivial solution for (p and x limit, Je and Jo — > 0, 

besides the trivial solution: = and x = 0- This is the inversion method. Thus 
we can get the non-perturbative solution 7^ or x 7^ 0, if it exists at all, by the 
inversion method ||4^ 



5.2 case (I) 



In the case of Ji > 0, J2 < (Jg > Jo), eq. ( [4.32| ) shows that, up to the leading 



order of g = l/Nf, and x are written as functions of two sources Je, Jo- 

:= Kl-iPJ! + QJ^) + 0{j'), 

X ■-= LJq - R{2jJo) + 0{J^), (5.5) 

^ Here we have adopted the UV cutoff Af to make the order parameter dimensionless. How- 
ever this is not a unique choice. We can take other quantities which play the same role, e.g., 
a := Nfe^. The result on the existence of the phase transition is unchanged irrespective of what 
quantity we may choose to define the dimensionless order parameter, as already pointed out in 
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where we have defined the dimensionless source: 



Je,o := (5.6) 

Here we find that the cross term JeJo does not appear in cf) and that , terms 
do not appear in x- In what follows O(J^) implies 0{Jl), 0{JeJo), 0{J^Jo) 
or 0{Jg). The coefficients K, L, P,Q, R are written as power series of g, e.g., 
K = H^=Qg^Kn. In this paper we pay attention to the first two terms K^^Ki 
up to 0{g) and similarly for other coefficients, L,P,Q, R. Our definition of 0, x 
leads to 

Ko = Lo = 1, Po = Qo = Ro = |, (5.7) 

as already shown in section 3. 

Now we invert the above set of equations with respect to two sources, Je and 
Jo- We assume the inverted form as 

I = T,<f) + A<p^ + Bx' + C<j)x + 0{^'), 

Jo = ToX + Dcf' + Ex^ + Fclyx + Oi^^), (5.8) 

where 0{^^) denotes C(0^), C(0^x), 0{(t)x^) or 0{x^). 

From the condition for the inverted series eq. ( |5.(j| ) to be consistent with the 
above equation eq. ( |5.5| ), the coefficient is determined. Hence the inverted series 
in the case (I) is obtained as 

Jo = ToX + F(l)x + 0{v^), (5.9) 

where 

Te = K-\ 

To = L-\ 

A = PtIK-\ 

B = QtIK-^ = QK-^L-\ 

F = 2RTeToL-^ = 2RK~^L-\ (5.10) 

together with C = D = E = 0. Therefore, up to the leading order of g, we obtain 

Te = l-K,g + 0{g'), 
To = l-L^g + 0{g^), 

A = Po + (Pi-3iriPo)^? + 0(/), 

B = Qo + iQi-K^Qo-2L,Qo)g + 0{g^), 

F = 2Ro + 2{R,-K,Ro-2LMg + Oig^). (5.11) 
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5.3 case (II) 

Next we consider the case where Ji, J2 > (Je < Jq)- In this case eq. ( [4.37| ) leads 
to 

:= KJ,-U{2jJo) + 0{J^), 

X := LJo-{SJl + TJl) + 0{J^). (5.12) 
Using the consistency condition, we obtain the inverted equation: 

1 = ToX + D(l)^ + Ex' + Oiip''), (5.13) 

where 

Te = K-\ 
To = L-\ 

C = 2UTeT,K-\ 
D = SrlL~\ 

E = TtIL-\ (5.14) 

and A = B = F. Since 

Ko = Lo = 1, So = To = Uo = y, (5.15) 

we obtain 

Co = 2Do = 2Eq. (5.16) 



6 Chiral symmetry breaking in pure QED3 

In this section we consider the pure QED3 without the bare Chern-Simons term. 
The induced Chern-Simons term H^q {k) is also neglected in this section as well 
as the bare Chern-Simons term. The Chern-Simons terms are included in the 
subsequent sections based on the result of section 4. This section is an independent 
one which illustrates the inversion method in a slightly simplified setting. 



6.1 single order parameter in QED3 

As a special case of the inversion scheme presented in the previous section, we 
introduce the single source term £j = Je^{x)'^{x) and consider a single order 
parameter which has the power-series expansion with respect to the expansion 
parameter g {g = or l/Nf): 

00 

= h[Je] = J2g^MJe). (6.1) 

n=0 
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After inversion, we will get 

oo 

Je = m = E a'^fM)- (6.2) 



m=0 

The coefficient function fmi4>) is related to (pn{Je) from the consistency condition: 

= Hfm. (6.3) 

In QED3 without the Chern-Simons term [(3 = 1,6 = 0), it is shown in the 
following that (p obeys 

:= KJ,- PJl + 0{Jl). (6.4) 
The inverted equation is obtained as 

Je = r0 + A02 + O(03), (6.5) 
where the consistency condition leads to 

r = K-\ 

A = Pt^K-\ (6.6) 

Hence, to the leading order of g, 

r = l-K,g + 0{g^), 

A = Po + {Pi-3K,Po)g + Oig''). (6.7) 

Even in the limit Je 0, (p in eq. ( |6.5| ) has an non-trivial solution (p = —tA^^ > 
when r < (for A > 0). Moreover, the effective potential is obtained as 

^(0) = + f 0', (6.8) 



from the relation 

SV{<P) 



Je. (6.9) 



Hence the non-trivial solution > is energetically more favorable than the 
trivial one = when r < (as long as A > 0). Therefore the chiral symmetry 
is spontaneously broken for g > gc = where gc is identified with the critical 
point. 



6.2 quenched QED3 

In this subsection we consider the chiral symmetry breaking in the quenched 
QED3. Here the term "quenched" implies that all the fermion loop correction 
to the photon propagator are neglected. This situation is realized by taking the 
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limit Nf — i> in an appropriate way. Therefore we estimate the critical point in 
the lowest order of the coupling constant e instead of 1/iV/. As shown in section 
2, we need to calculate 



1 r (fk „,r,\,,. d 
2J (2^ 



1 f^v ,2 (9 
dk- 



47r2 7o P dJe 
To the lowest order O(e^), we obtain 



n^.'^^)]. (6.10) 



1 r d^k d 

-I 

e 



[ ^/)(o)(A;)^n«(A;) 

2 7 (27r)3 ^""^ 'dJe ^'^^ ' 



where we have temporally introduced the infrared (IR) cutoff besides the ultra- 
violet (UV) cutoff Ap. For a while, we put (3 = 1. In the quenched limit, all the 
loop corrections are neglected and hence the photon propagator reduces to the 
bare one and the order parameter is defined by taking the limit Nf — > 0. Then 
we have 

^2 /iTrQiTraX / ^2 



lim —^^^ = jj 1 + ^ In ^ ) - PJ2 + 0( J^). (6.12) 



After inversion up to C(e^), we get 

J, := ^ = f 1 - ^ In A<P' + 0{<P'). (6.13) 

Even when Je = 0, there may exist a non-trivial solution 7^ besides the trivial 
one, = 0: 

= f^ln^-lV (6.14) 

Actually, defining the dimensionless inverse coupling constant (3e ■= A//e^, we see 
there is a non-trivial solution when /3e < /?e for a critical value: 

/?::=^ln^. (6.15) 
8 /i 

This result however shows that (3^ ^ 00 as fi I 0. Naively this may be interpreted 
as follows. Q In the quenched QED3, the chiral symmetry is spontaneously broken 

^ Note that we obtain 

7:^Aa8 f^- 1 

which diverges as | 0. 
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for arbitrary values of the coupling constant 7^ 0. This agrees with the results 
of the Schwinger- Dyson equation approach and the Monte Carlo simulation 
of the lattice non-compact QED3 PB| . 



6.3 unquenched QED3 in the l/Nj expansion 

Next we consider the problem of chiral symmetry breaking in unquenched QED3 
based on the scheme of the 1/iV/ expansion. To the leading order in expan- 
sion, we wish to calculate 



27r2 7o 'dJe 

' k'dk-——^^^^. (6.16) 



2vr2 7o [Pk^-U?{k)] dJ, 

Up to 0{J^), this reduces to 



^ '''''k'dk--^-—8a(^--^]+0{J!). (6.17) 



271^ Jo [(3k'^ + ak] \k vr P 

Therefore we obtain 

After inversion, we get 

l=h-jJ-U + Act? + 0{ct>'), (6.19) 

withQ 

m^IU^^'-^ r dk-^ = lu(l + (6.20) 

In the leading order of the 1/A^/ expansion, we have shown the existence and 
finiteness of the critical number of flavors N'j, below which the chiral symmetry is 

^ Note that we obtain 

Pi - / dk- 



2 7r2 TT k{pk + a) ' 

Naively this diverges as /i | 0. The infrared divergence in Pi can not be removed in this 
method. This term may induce the fluctuation-induced first order transition as discussed in the 
final section. 
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spontaneously broken and above which the chiral symmetry restores. This resuh 
should be compared with the previous analysis in the quenched case. The reason 
why we can obtain a finite critical number is that the expansion alters the 

infrared behavior of the gauge boson propagator from to l/{ak) due to the 
existence of Ii^^\k) in the denominator. This eliminates the infrared divergence 
encountered in the previous subsection. 

Moreover, even in the Chern-Simons limit /? — >■ 0, i.e., without the kinetic term 
for the gauge boson, the finite critical number of fiavors is obtained: N'j{f3 = 0) = 
2. When /3 7^ 0, the naive "continuum" (infinite cutoff) limit A/a ^ 00 leads 
to A^j — > 0, i.e., no breaking of the chiral symmetry. However the first viewpoint 
explained in the introduction that the dimensional quantity a gives a natural scale 
of the theory implies that the relevant physics comes essentially from the region 
/c ~ a and the region k ^ a does not give the essential contribution. Hence we 
should normalize the physical quantity by a (instead of Aj) and set the upper 
bound of integration as Ap ~ a. This standpoint gives a finite and nonzero value 
for N'j. 

7 Effective potential and stability 

7.1 chiral- symmetry violation (case (I)) 

We define the dimensionless effective potential r[0, x] by 

f u 1 I- 
r[0,x] = ^^^, (T.i) 

in conformity with the definition of the dimensionless order parameter (|5.1|) . We 
consider the effective potential 

f [0, X] = ln<P' + \roX' + \Act>' + \fcI^x' + 0{v'), (7.2) 
which satisfies the following relation associated with the case (I): 

= Je = r.0 + V + ^i^X^ + O(^^), 

= Jo = T.X + F<t^X + 0{^'). (7.3) 

ox 

It is easy to see that, even when Jg = = 0, a pair of eq. ( [7.3|) has non-trivial 
solutions besides the trivial one = x = 0. First of all, we observe that a pair of 
solutions such that x 7^ and = does not exist, which is consistent with the 
assumption Jg > Jo in case (I). 

There are two types of non-trivial solutions. One type of non-trivial solutions 
is given by 

x(i) = 0, 0(^) = -TeA-i > 0, (7.4) 
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where Xe < and A > 0. At this solution, the effective potential takes the value: 

r[cf>('\x^'^] = lr!A-'<0, (7.5) 

which is lower than the trivial one: r[0 = O,x = O] = O. This solution implies the 
spontaneous breaking of chiral symmetry and no spontaneous breaking of parity. 
Another non-trivial solution is given by 



= -roF-\ x''' = (7.6) 

The solution 0*^^) is positive for Tq < and x^^"* is real and positive if Tq < and 
FT(, < Ato, as long as A, F > 0. This solution implies that the chiral symmetry 
and the parity are simultaneously broken. It gives the effective potential 

n¥'\ X^'^] = \tIF-\?,Ft, - 2At,) < ]:t!AF-^ < 0. (7.7) 
o o 

Then this solution also gives lower effective potential than r[0, 0] = 0. 

Which solution gives the absolute minimum? If the second solution is the 
absolute minimum, the chiral symmetry and the parity is spontaneously broken 
simultaneously. However we can show that the first solution is energetically fa- 
vorable: 

f[0«,x«]<f[</)(^),x(')]<O, (7.8) 

since 

f[0«,x«]-f[0(2),x(^)] 

= [2|Ar,|3 + \Ft,\' - 3\Aro\''\Fn\]/[6A^F'] < 0, (7.9) 

in the allowed region of parameters. A, F, Tg, Tq. 

Now we consider the stability of the solution around the respective stationary 
point: (j) = X = X^"'' which |^ = and |^ = 0. The effective potential can 
be expanded around the stationary point as 

r[0,x] = n^^^\x^-^] + 5x)Hr[^^-\x^'^^] (^^) + OiiS^r), (7.10) 

where 6(p denotes deviation 6(j) or 6x from the respective stationary point: 50 := 
— 0("\ Sx ■= X ^ X^"'' ^iid -ffr[0, x] is the Hessian matrix defined by 

^r[0,x]:= '^]. (7.11) 

For the above effective potential ( |7.2| ), the Hessian reads 
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Around the first solution, 



Note that —t^ > and {Ato — Ft^)/A > 0. Hence trifp > and the Hessian is 
positive definite, i.e. det i^r > 0. Then the first solution x^^'^ gives the local 
minimum. 

Around the second solution, 

^,|^<^x'=)]^((^--^,t°>''^ ). (7.14) 

The Hessian is negative definite, since det Hj- = —{Fx^'^^Y < 0- Hence the second 
solution (f)'^'^\ x^^) corresponds to the saddle point. 

Thus the first solution give the absolute minimum. Therefore we can conclude 
that there are no stable (absolute) minima at which both (p and x ^"^^ simultane- 
ously non-zero. This conclusion agrees with the result of Semenoff and Wijeward- 
hana [^. Moreover we have shown up to 0{1/Nf) that the chiral symmetry is 
broken spontaneously for Nf < NJ = Ki, while the spontaneous breaking of par- 
ity does not occur. In the broken phase the non-trivial order parameter behaves 
near the critical point as 

0(1) = -TeA-i = {K,/Nf - 1)P,-' + 0{1/Nf), (7.15) 

where the value Ki is given by eq. ( [4.351) . Note that Ki is always positive and 
monotonically decreasing in (3 A/ a and has a maximum value 2 at the Chern- 
Simons limit: (3A/a = 0. Hence N'^ has a finite upper bound 2. For example, 
when (3A/a = 1, Ki = 21n(2) = 1.4. 

7.2 parity violation (case (II)) 

Next we consider the case where Ji > 0, J2 > (Jg < Jo)- This case corresponds 
to the effective potential: 

f[0, x] = \rect^' + \rox' + \{Ex' + \cxct>') + 0{^^), (7.16) 
since this gives the following relation for the case (II): 



5(p 

- Jo = roX + ^<P' + Ex' + 0{v'). (7.17) 

One non-trivial solution is given by 

0(1) = 0, x^'^ = -ToE-^ > 0, (7.18) 
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where Tq < and E > 0. This solution gives the effective potential 

r[0«,X«] = ^r^E-2<O, (7.19) 

which is lower than r[0 = 0,x = 0] = 0. This solution implies the spontaneous 
breaking of parity and no spontaneous breaking of chiral symmetry. 
Another non-trivial solution is given by 



X''^ = -reC-\ = ^1^^^^^^. (7.20) 

The solution x*^^) is positive for Tg < and x^'^^ is real and positive if Xe < 0, Ctq < 
ETe, as long as C, > 0. This solution gives the effective potential 

r[0('\ X^'^] = ]:r',C~'{3Cr, - 2Et,) < 1:t!C-'E < 0, (7.21) 

D 

which is lower than r[0, 0] = 0. 

We can show that the first solution is energetically favorable: 

f[0«,x(^^]<f[0(^),x(')]<O, (7.22) 

in the same way as the previous case (I). 

For the above effective potential (|7.16| ), the Hessian reads 

HrlM-r-^^"^ ^J). (7.23) 



C4> 2E 



Around the second solution, 



The Hessian is negative definite, since det Hr = — (C^*^^-*)^ < 0. Hence the second 
solution 0*^^^ x^"^^ corresponds to the saddle point. 
Around the first solution, 

= ^l). (7.25) 

Note that E > 0, {Et^ — Cto)/E > 0. Hence triYp > and the Hessian is positive 
definite, i.e. det i^r > 0. Then the first solution 0^^^ x^^^ gives the absolute 
minimum. Therefore we can conclude that there are no stable (absolute) minima 
at which both (p and x ^^re simultaneously non-zero. 

Hence we are tempted to conclude that up to 0{1/Nf) the parity symmetry 
may be broken spontaneously for Nf < N'j = Li, while the spontaneous breaking 
of chiral symmetry does not occur. If so, the parity-odd order parameter should 
behave as 

^(1) = -roE-' = {Lr/Nf - l)S^^ + Oil/N}) > 0. (7.26) 
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However note that the value Li given by eq. ( [4.39| ) has a maximum value 0.1208 
(at X := (3 A/ a = 13.3) which is extremely small. Above which Li decreases 
monotonically and goes to zero at x — oo, while Li becomes negative below x = 
3.92 and reaches to Li = —2 at the Chern-Simons limit: x = 0. This seems to show 
that the spontaneous breaking of parity does not occur actually. The absence of 
spontaneous breakdown of parity in (2+l)-dimensional gauge theories is consistent 

j, Eol, O, 02], p6i and general consideration by 



with results of other works 38 



Vafa and Witten g|. 



8 First order transition 



8.1 effective potential in the presence of bare Chern-Simons 
term 

In this section we study the effect of the bare Chern-Simons term. In what follows 
we restrict our consideration to small 6 for simplicity of calculation. It does not 
cause any problem to consider the large 9 case, in principle besides the calcula- 
tional complexity. Taking into account eq. ( [4.31| ), we obtain from eq. ( [4.26 ): 



X 



aeJo + KJ, - {PJl + QJl) + 0{ej\ J^), 
be + aeJe + LJ, - R{2jJo) + 0(^J^ J^), 



where ae and be are quantities of order 9/Nf given by 

TT^ e 32 r^v d^k 



ae 



2K) Nf TT 



a 



2A} Nf 



(27r)3 (/?A;2 + aky 
d?k Aak 



O 



{2T^f {(3k^ + akf \Nf 




(8.2) 



The inverted equation is given by 



Je 
Jo 



Deep + T,<P + V + Bx^ + C( V, <^') 
Ce + Eect) + r,x + Fcj^x + 0( V, <^'). 



(8.3) 



From the consistency condition, the relations among the coefficients are obtained 
in the same way as before. For example, we obtain: 



e 
To 



Q-'ae = -L-% ~ 0{e), 
aeR-^To = aeK-^L-\l - a^R-^L-^y^ 
aeL-We = aeK-^L-\l - ajR-^L-^)-^ -- 
K-\l - alK-^L'^)'^ = K'^ + 0{e^), 
L-\l - alK-'L-^)-^ = L-^ + 0{e^), 



-- aeK-^L-^ + 0(6^) 



•4) 



where we have used R = Q = —[si — S2 — (ti ~ ^2)], see eq. ( [4.26| ). In the limit 



9 



0, we see Ce, De, Eg 



o{e) 



0, and Tp 



K 



-1 



L , as expected 
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from the result of section 5. As in section 7, the effective potential is obtained by 
integrating eq. 



m X] = Cex + Ee<Px + ^^e0' + \r,X^ + \acP^ + \fcPx^ + 0{eip\ if'), (8.5) 
so that the effective potential satisfies 

5f ~ 5f ~ , , 

— = Je, — = Jo, 8.6 
b<p bx 

where we have used Eq = Dg. 

Let 0*, X* denote the stationary point of the effective potential r[0, x], i-e., 



Sr 6T 

=</> :X=X ox 



,x=x 



0. (8.7) 



In order to specify the location of the stationary point, we put r := x*/4'*- From 
the first equation of eq. (P?B|), we get 



^ A + Br^' 

On the other hand, the second equation of eq. (|8.3| ) leads to 



Ee + Tor ± ^{Ee + Tpr f - AFCer 



2Fr 

Equating two values for (p* in the above two equations, we get the algebraic 
equation for r which is 4th order in r. Since we restrict our attention to small 6, 
r is obtained up to 0{6) as 

rs - A^^^ . an (8.10) 

since Ce.Ee^ 0{9). 

At the stationary point we define an order parameter (f by 

^ X* = 7^-^V. (8.11) 



^ (l+r2)i/2^' ^ (l + r2)i/2 

In the case of = 0, = and hence (0*, x*) = if, 0) which is consistent with 
the result of section 7.1. In the plane x* = ^o4>* including the stationary point 
{(f)*,X*), the effective potential behaves as (see Fig. 4) 

This is a generalization of eq. ( [7.2| ) to the case of ^ 7^ 0. 
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8.2 1st order transition 



In the presence of the bare Chern-Simons term, we are forced to investigate the 
behavior of the effective potential of the form: 

= A^{Nf, e)ip + \A^{Nf, e)ip^ + ^^(iv^, e)ip^ + o{eip\ y.^), (s.is) 

where all the coefficient Ai{i = 1,2, ...) are functions of the bare parameters: Nf 
and 9. 

First we consider the situation in which the Chern-Simons term is absent: 
Ai = and ^3 > 0. For positive ^2 > (point A), the effective potential is 
a convex function of (p in the region if > and has an absolute minimum at 
!f = 0, see Fig. 5 For negative ^2 < (point B), it has a minimum at (p = '■ = 
—A2/AS > 0. This region of parameters corresponds to the chiral symmetry- 
breaking phase. As the parameter A2 passes through A2 = from the positive 
side to negative one, the absolute minimum at ^9 = moves continuously to the 
absolute minimum at ip = 7^ 0. Therefore this theory exhibits a second order 
phase transition at A2{Nf = NJ, 0) = 0. This determines Nj, the critical number 
of flavors, below which {Nf < NJ) the chiral symmetry is spontaneously broken. 

Now we can show that the non-vanishing Ai changes the type of the phase 
transition drastically. In this case, there occurs a first order phase transition at a 
certain value of Ai. In order to show this, we start from a point B in the broken- 
symmetry phase of the phase diagram {Nf,6) such that Nf < Nf and ^ = 0. 
At this point, Ai{Nf,e = 0) = and A2{Nf,e = 0) < 0. As the Chern-Simons 
coefficient 6 grows from zero, the shape of the effective potential T[ip] changes 
as shown in Fig. 5. For small 6 (point C), Ai > is also small and T[ip] has 
an absolute minimum at = ipm{Nf,6) away from ip = 0. As Ai is increased, 
the value r[(y9m] at = ipm increases monotonically and approaches to zero at 
Ai = AI (point D). Therefore there is a certain value of Ai such that T[ipc] = at 
Ai = AI and above which {Ai > AD the absolute minimum moves to the origin 
discontinuously (point E). This is nothing but the first order phase transition 
point. At the critical point, the effective potential takes the form: 

= ^AsiNf, OMip - v..)^ + 0{v'). (8.14) 

Therefore there is a critical value 6c{> 0) of Chern-Simons coefficient at which 
the discontinuous first order phase transition occurs in the region Nf < Nf. By 
comparing eq. ( |8.14| ) with eq. ( |8.13D , the critical value of 9c = Oc{Nf) as a function 



of Nf is determined: 9c satisfies the following equation: 



and (pc is given by 



3 A2iNf,9c) _ 3A,iNf,9c 
AA,iNf,9c) \ A,iNf,9c: 



:> 0). (8.16) 
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Here note that rg is proportional to the Chern-Simons coefficient and hence 
6 = imphes Ai = 0. Therefore the ffist order phase transition in Maxwell- Chern- 
Simons theory with fermions occurs only when there is a bare Chern-Simons term, 

Therefore, comparing eq. ( ^.12| ) with eq. ( ^.13| ), we find that 



A2 

As 



(8.17) 



and hence the critical line for the first order chiral phase transition is given by 



N. 



f , 



in the neighborhood of the critical point {Nf,6) = {N'^,0). This result should be 
compared with that of Schwinger- Dyson equation approach 



exp 



TT 



IA9) 



The result eq. ( |8.18| ) is reasonable, since it is obtained up to the leading order of 
1/Nf expansion. To the lowest order of 1/Nf expansion, we can not recover the 



essential singularity behavior even if it is correct, see |49 . 



9 Conclusion and discussion 

In the leading order of expansion, we have constructed the gauge-covariant 
(gauge-parameter-independent) effective potential in terms of two order parame- 
ters, {^'^) and {^T"^), in the Maxwell and/or Chern-Simons theory coupled to 
fermionic matter. From the viewpoint of the stability of the solution around the 
stationary point of the effective potential, we have shown that the spontaneous 
breaking of both the chiral-symmetry and the parity can not occur simultane- 
ously. This result in the absence of the bare Chern-Simons term agrees with that 
obtained only in the Landau gauge by Semenoff and Wijewardhana In (2+1)- 
dimensional gauge theories without the bare Chern-Simons term {6 = 0), there 
may exist a finite number of critical flavors such that for Nf < Nj the chiral 
symmetry is spontaneously broken. On the other hand, the spontaneous break- 
down of parity does not occur. This conclusion agrees with the previous analyses 
[P3, pO, HO, HTI, H^, EBI and general consideration P3[ . We wish to emphasize that 



all the results obtained in this paper are gauge-parameter independent. 
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Moreover, in sharp contrast to the Schwinger- Dyson equation approach of 



QED3 up to now (see e.g. 1|22|), we have included the effect coming from the 
vacuum polarization to the photon propagator exactly, up to the leading order of 
l/A^j expansion. In the Schwinger- Dyson equation approach for QED3 so far, the 
vacuum polarization effect to the photon propagator has been included through 
the expression Ilip\k) = ak which is obtained in the massless fermion limit of 
Il'^\k). In order to include the vacuum polarization effect from the massive 
fermion in the Schwinger- Dyson equation approach, we must solve the Schwinger- 
Dyson equation for the photon propagator simultaneously with the Schwinger- 
Dyson equation for the fermion propagator. In the four-dimensional case such 



a calculation has been performed in the strong-coupling phase of QED4 |52]. 
Therefore the result obtained in this paper based on the inversion method goes 
beyond the Schwinger- Dyson equation approach for QED3 done so far in a sense 
just mentioned. However our analysis shows that the inclusion of the massive 
fermion does not change the critical number of flavors, if any. 

The critical value Nj obtained here seems rather small compared with the 
result of the Schwinger-Dyson equation Nj = 3 ~ 4 and the Monte Carlo simu- 
lation, Nj = 3.5 ± 0.5. This will be due to the fact that the above calculation 
is restricted to the leading order of l/Nj expansion, although the calculation is 
systematic. This is also the case for QED4, see P9|. To estimate more accurately 



the critical flavor, we need to perform the higher order calculation. 

The order of the chiral phase transition is of the second order in the absence of 
the Chern-Simons term. This is consistent with the result of However, in the 
presence of the Chern-Simons term, we have shown that the chiral phase transition 
turns into the first order transition, as discovered firstly in the Schwinger-Dyson 
equation 



Now we come to the stage of discussing the meaning of the "continuum" limit 
(of removing the cutoffs) in the non-perturbative sense The point (iVj, 0) 

is the only one point which exhibits the continuous second order chiral transition. 
Therefore, in order to obtain the continuum theory with a finite fermion mass, 
the bare parameters must be adjusted so as to approach the critical point {NJ, 0). 

It should be remarked that the bare Chern-Simons term plays completely 
different role from the induced Chern-Simons term. It is the bare Chern-Simons 
term which causes the first order transition, while the induced Chern-Simons term 
does not leads to the first order transition, even if the Chern-Simons term may 
be induced by radiative corrections at one-loop. The Chern-Simons coefficient 6 
is not subject to renormalization according to the Coleman-Hill theorem and 
the beta function of the renormalization group is identically zero: f3{6) = 0. The 
continuum theory obtained from the massive phase in this continuum limit may 
have no Chern-Simons term in the renormalized sense. 

In this paper we have neglected the possibility of the fluctuation-induced first 
order transition. If this happens in this theory, the point {NJ, 0) will be no longer 
of the second order transition point. To clarify this problem, we need to perform 
more careful analysis on the infrared behavior of the theory. However the absence 
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of such a first order transition has been reported recently in 

Finally we remark the possible effect coming from higher orders in ex- 
pansion [^]. Our analysis of the effective potential up to the leading order of 
expansion shows 



(9.1) 



in QED3 without the Chern-Simons term. In the Schwinger-Dyson equation ap- 
proach the chiral order parameter as well as the dynamically generated 
fermion mass exhibits the essential singularity at the critical point iVj = A^^: 



3/2 



a 



2 



\ (X ) 



exp 



(9.2) 



In the four-dimensional case, we have shown how the result obtained in the 
Schwinger-Dyson equation approach for QED4 can be reproduced from the higher 
order calculation in the inversion method. For this, it is helpful to enlarge the 
model so as to include the four-fermion interaction. We can expect that this 
scenario also holds in the three-dimensional case, see e.g., However, even if 



this conjecture is correct, we will need the infinite order result of l/iVj expansion 
to recover the essential singularity exactly. A similar feature can be seen in the 
expression of the critical line between eq. ( p.l^ ) and eq. (|8.19|) . The result A^j = 2 
up to the leading order for the Chern-Simons theory [fi = 0) with Nf flavors of 
4-component fermions is not so large that the reliability of 1/Nf expansion may be 
not very definite. To see the improvement of the expansion by higher orders, we 
need to calculate the next-to-leading order. This is the subject of the subsequent 
work. 

One more comment is on the choice of order parameter. In the analysis of this 
paper, we did not consider the order parameter 

a; := ^ 1 d'x{e^'''A,{x)d,Ap{x)). (9.3) 

This order parameter u as well as x is able to signal the spontaneous breakdown 
of parity after taking the limit 6^0. In this case we must obtain the effective 
potential as a function of three order parameters 0, x a^id u. This is possible 
in principle, but is somewhat tedious. To discuss the connection of our result 
with the anyonic superconductivity, we must include the chemical potential into 
this theory. This will be done in the forthcoming paper. The first order phase 
transition induced by the topological term may have other important implications, 
see \^%. 
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A Clifford algebra and its representation 

In the Minkowski space-time, the gamma matrices satisfy the Chfford algebra: 

rMlM + lWM = '^9'''. (A.l) 

which is satisfied for example by 
with Pauli matrices: 

In the euclidean space, we define the gamma matrices satisfying 

7^7^ + 7^7^ = -25'^^ (A.4) 
by 1e ■= -ilM,lE ■= Im^Ie ■= 7if, i-e., 



Note that 



and 



(7a' = -(0 5) (f = 0,l,2) (A.6) 



Then we define 



7|:=(° n)> (A.8) 



7|:=7l7kl;7l=(? V)' ^^'^^ 



r-lh%-(i \]. (A.10) 



and 

.0 -/ 
For the projection operator: 
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we obtain the following trace formulae: 



and 



tr[x±] = 


= 2, 






tr[7p7^7pY±] = 

L / I2j I Hi I xZ//v— J 




tr[7K7K7K7KX±] = 


= 2( ^./i/^rv/S 


tr[r] = 


= 0, 


tr[7^7^T] = 


= T25j^^, 


tr[rx±] = 


= ±2, 


tr[7£7£X±] = 


= -25;,,/, 


tr[7^7E7EX±] = 


~ ~F2c^,yp, 


tr[7^7^7S7^X±] = 





(A.12) 



(A.13) 



B Generating functional 

For the bosonic field, A, the functional integration is performed as 

/ VAexp[--AD-^A] = (DetD^^)-^/^ = exp -InDetD . 
t/ 2 _ 2 _ 

Since the full boson propagator D satisfies the Schwinger-Dyson equation: 



we have 



D = - n]-i = [1 - L>(°)n]-iL>w, 



lnDet[D] = Trln[D] 

= Tr ln[D(°)] - Tr In[l - D^^^U] 

CO 1 

= Trln[L>(°)] + ^-Tr[(i^("'n)^ 



n=l 



n 



Then the derivative takes the form: 



8 r9 °° 1 

oo r a 



n=0 

Tr 
Tr 



'i_^(o)n)-iD(o)An 



(B.l) 



(B.2) 



(B.3) 



(B.4) 
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Especially, for the gauge field, we obtain 

d 



InDctiD] 



(B.5) 



On the other hand, for the fermionic field, we have 

J Vi^Vi^ expfV'^" V] = Det^-^ = exp [- In Det5] . (B.6) 
Since the full fermion propagator S satisfies the Schwinger-Dyson equation: 



S = [^(°)" - E]-i = [1 - ^(o)E]-i^(o), 

we obtain in the similar way to the bosonic case 

lnDet[5] = Trln[S] 

= Trln[5(°)] -Trln[l-5(°)E] 

OO 1 

= Trln[,5(°)] + ^-Tr[(5(°)E)' 



(B.7) 



n=l 



n 



and 



d 



InDetf 



d 



Trln[5(°)] +Tr 



dJ ' ' dJ 

In the scheme of 1/iV/ expansion, we find 

e'^Nf = a. 



4-' 



(B.8) 



(B.9) 



a 



~ e 



N 



f 



(B.IO) 



Therefore, in the leading order of 1/Nf expansion, we have only to calculate 



B.ll) 



C Vacuum polarization 

C.l vacuum polarization tensor 

The vacuum polarization tensor at one-loop is obtained from 

1 1 



nW±(^;-r„) = -e 



' I ^tr 



7m- 



-Iv- 



-x± 



— — e 



-\- )k — m p — m' 
2 1 f^^P tr[7^(j^+ /; + TO)7^(j6 + m)x±] 



J (27r 



(27r)3 [(p + A;)2 + m2](p2 + ^2) 



(C.l) 
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Hence the vacuum polarization tensor is decomposed into two parts: 

d^p tr[7^(/>+ /;)7,.i)X±] + mhT[-f^'y^x±\ 



and 



—e^m 



(27r)3 (p2+^2)[(p+^)2+^2] 

d^p tr[7^(i>+ /;)7^x±] +tr[7^7^, 



(27r)3 (p2 + ^2)[(p + ^)2 + ^2] 

The parity even part is calculated as 

d^p tr[7^(i)+ /c)'y^^x±] + mhT[-f^-f^x±\ 



(C.2) 



(C.3) 



— e 



(27r) 



(p^ + m^)[(p + ky + rn?] 



kij,ki/\ e 



2 ^1 



- ( - ^ 1 ^ 



x(l — x) 



kiikj. 



V '^'^ F ; 27r 
On the other hand, the odd part reads 



1 k^ - 4m^ k 
- vm^ H ; arctan — ^= 

2 4A; 



(C.4) 



/ 



= ^2e^me^ypkp dx J 



d^P tr[7^(i)+ /:)7i^X±] +tr[7^7^i>X±] 

(27r)3 (p2+^2)[(p + ^)2 + ^2] 

d'^p -f'^^ixupkp 

(27r)3(p2+^2)[(p + ^)2 + ^2] 



4^ e^i^P^P 



(27r)3 [p2 -|- 2x/c • p + xk^ + rn?]'^ 
dx 



— — ^ui/pkp — arctan 
47r k 



[m'^ + x{l — x)k'^Y/'^ 
k 



2vm2 



(C.5) 



where we have used tr[7^7i,7pX-|-] = =F2e/ii/p and the Feynman parameter formula 
AB — Jo "-^ [ 



AB JO "^-^ [xA+(l-x)B]^ ■ 

For the 2-component fermion, we have 



2 1 J I H ; arctan 



k 



2 J 



6^ k 
nff{k;J) = — Jarctan^, 



and the the derivative with respect to the source J is given by 



_d_ 
dJ 



e 

An 



sgn{J) 



8J2 



fc2 + 4 J2 k 



4J 

— arctan 



2|J| 



4|J| 



A;2 + 4J2 A; 



2 A; 
+ T arctan 



2|J|, 



(C.6) 
(C.7) 

(C.8) 
(C.9) 



where k :— Vk^. 
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C.2 power-series expansion (I) 

Using the mathematical identities: 



X 



arcsm 



— arctanx = ^ 



{ 2n-l 



X 



2n-l 



TT 



- arctan(a;-^), (C.IO) 



we obtain the following expansions. 
For < 1, we get 



16 



TT 



n=l 



- 1 ^ A; 



4\ H 



(C.ll) 



and 



-J 



-J 



1- sgniJ)- 



n-l 



TT 



n=l 



2n 



_/^\2n-l 
V k ' 



TT /C 



9J 



n5!^(A;; J) 



e 
2" 



r k 



1 2 °° nil 



2n 



(C.12) 



(CIS) 



|,ng)(M) 



e 

T 
4 



4 oo 

l-55n(J)-^(-ir-^ 



n 



TT 



n=l 



2n 



_/^N2n-l 



i_ + + o 



(C.14) 



C.3 power-series expansion (II) 



For 1^1 < 1, we obtain 



I 2JI 



,2 oo 



Il'^\k-J) = sgn{J)-kJ2{-ir 

n=l 

-.2 



A: 



27r ±1' ' An^-V2J 



sqn(J) — k 



1 k 

-37 + ° 



3 



,2n-l 



(C.15) 



^2 oo 



n-l 



27r 



rt=l 



s^n(J)— A; 



2n- 1'2J 

2 



(C.16) 
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and 



dJ 



^2 oo 



n 



\2n 



TT 



n=l 



2n 



sgniJ)- 



37^ + ^ 7 



r2j 

\ 4 



(C.17) 



9J 



^2 oo 



2n-l 



n=2 



(C.18) 
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Figure Captions 

Figure 1: Schematic phase diagram {Nf, 9). The sohd hne denotes the phase transition 
hne below which the chiral symmetry is spontaneously broken. 



Typical points in the phase diagram 
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1st order transition 
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chiral symm., parity B 



Figure 2: Dynamical fermion masses as functions of 9 for a pair of 2-component 
fermions. 



Figure 3: Feynman diagrams needed to evaluate the effective potential up to the lead- 
ing order oi 1/Nf expansion. The sohd line corresponds the fermion propa- 
gator and the broken line to the gauge boson propagator. 



Figure 4: The effective potential in terms of two order parameters (f) = (p, and x = V^o 
in the chiral-symmetry breaking phase. When ^ = (the point B), the 
effective potential has a minimum at (92*, 0). When ^ 7^ 0, the location of 
the stationary point moves to (</?*, </?*) where (fil and </?* 7^ 0. 



Figure 5: The shape of the effective potential corresponding to the respective point, 
A, B, C, D, E, in the phase diagram of Figure 1. The phase transition from 
AtoB corresponds to the second order one, while the phase transition at 
D is the first order one. 
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